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In recent years, the Ferm/-LAT gamma-ray telescope has detected a population of over 160 
gamma-ray pulsars, which has enabled the detailed study of electromagnetic radiation from pul¬ 
sars at energies above 100 MeV. Further, since the surprising detection of the Crab pulsar in very 
high-energy (VHE; E > lOOGeV) gamma rays by the MAGIC and VERITAS collaborations, 
there has been an ongoing effort in the astrophysics community to detect new pulsars in the VHE 
band. However, the Crab remains the only pulsar so far detected in VHE gamma rays, raising 
the question of whether or not the Crab is unique and also making it more difficult to constrain 
model predictions that attempt to explain the emission. Presented here are recent VERITAS re¬ 
sults from observational campaigns on the brightest northern-hemisphere high-energy gamma-ray 
pulsar Geminga and the missing link binary pulsar PSR J1023H-0038, which have both resulted 
in upper limits on a possible VHE flux. These limits are placed into context with the current 
theoretical framework attempting to explain the origin of the gamma-ray emission from pulsars. 
Additionally, future plans for pulsar observations with VERITAS will be briefly discussed. 


The 34th International Cosmic Ray Conference, 
30 July- 6 August, 2015 
The Hague, The Netherlands 


* Veritas, sao.arizona.edu 
' Speaker. 


© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








Geminga and PSR J1023+0038 results with VERITAS 


G. T. Richards for the VERITAS Collaboration 


1. Introduction 

The field of pulsar astronomy has entered a new era in the time since the launch of the Fermi 
Large Area Telescope (LAT), which has now succeeded in detecting over 160 pulsars in the high- 
energy (HE; E > 100 MeV) gamma-ray band* . These detections have greatly facilitated the detailed 
study of electromagnetic radiation from pulsars at the highest energies. The spectral energy distri¬ 
butions (SEDs) that have been seen for pulsars in gamma rays are all well-characterized by a broad 
curvature radiation component that originates due to the electrons and positrons that fill the mag¬ 
netosphere and follow curved trajectories as they are confined to magnetic field lines. A curvature 
radiation component has a natural end where the emission becomes radiation-reaction limited [1], 
resulting in a spectral break at a few GeV [2]. The spectral “tails” of pulsars in the GeV band are 
well-described by an exponential cutoff as expected for curvature radiation [2], though statistics are 
sparse above ~ 10 GeV due to the decreasing sensitivity of the Fermi-LKF at these high energies. 

In recent years, a few pulsars have been detected at energies considerably above the HE spec¬ 
tral break. The Vela pulsar, the brightest steady gamma-ray source seen by the Fermi-LKT [3], has 
been detected from the ground by H.E.S.S.^ above 30 GeV (Gajdus, M., these proceedings) [4] and 
by the Fermi-LKT above 50 GeV [5]. The H.E.S.S. detection of the Vela pulsar makes it the second 
pulsar to have been detected in gamma rays from the ground. The first, the Crab pulsar, has been 
detected in the VHE band above 100 GeV by MAGIC [6] and VERITAS [7] and has posed a chal¬ 
lenge in that the curvature radiation scenario is not adequate for a full explanation of the gamma-ray 
radiation; the combined Thn7j/-EAT and VERITAS SED favors a power-law fit above ~ 10 GeV [7]. 
Current models that have been proposed to explain the VHE emission from the Crab pulsar include 
inverse-Compton (IC) scattering scenarios in the outer magnetosphere [8, 9, 10] and beyond the 
light cylinder [11, 12], and synchrotron radiation from magnetic reconnection events [13]. Whether 
or not the VHE emission from the Crab pulsar is exceptional with respect to the rest of the gamma- 
ray pulsar population remains to be seen. 

In this paper we review recent results from the VERITAS collaboration on the Geminga pulsar 
and PSR J1023-I-0038. The paper is structured in the following way: in Sections 1.1 and 1.2 
we review the Geminga pulsar and PSR J1023-I-0038, respectively; in Section 2 we discuss the 
VERITAS observations, analyses, and main results; and in Section 3 we present a brief overview 
of future pulsar science plans with VERITAS. 

1.1 The Geminga Pulsar 

Eormally known as PSR J0633-I-1746, the Geminga pulsar is located at the relatively close 
distance of ~200pc [14, 15] and is the second-brightest steady HE gamma-ray source in the sky. 
The Geminga pulsar has a spin period of ~240ms and a spin-down luminosity of 3.2 x 10^4 erg 
s^' [16]. It is the first GeV-emitting pulsar detected with no known radio counterpart [17] and was, 
in fact, first discovered in gamma rays by the SAS-2 [18] and COS-B [19] satellites. Though radio 
searches have failed to reveal any hints of emission, the Geminga pulsar is detected in UV [20] 
and X-rays [21]. Similar to the Crab pulsar in the high-energy band, Geminga shows a double- 
peaked phaseogram separated by a “bridge” of emission above the background, with only one of 

*https://confluence.slac.stanford.edu/display/GLAMCOG/Public-l-List-l-of+LAT-Detected+Gamma-R ay+Pulsars/ 

2http://www.mpg.de/8287998/velar_pulsar/ 
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the peaks remaining dominant above about lOGeV. As has been seen for other gamma-ray pulsars, 
the spectrum of Geminga above lOOMeV is can be described by a power law with an exponential 
cut-off [2], characteristic of the curvature radiation processes thought to dominate in the magneto¬ 
sphere. However, [10] reported that the spectrum is better described by a simple power-law above 
the break energy, bringing into question that the spectrum is best described by the exponential 
cut-off expected in a curvature radiation scenario. 


1.2 PSR J1023-I-0038 


The system PSR J1023-I-0038 is an eclipsing binary system located at a distance of 1370 ± 
40 pc [22] containing a millisecond pulsar (MSP) with a rapid rotation period of 1.69 ms orbiting 
a non-degenerate companion star every ~5 hr', found in a survey in 2008 by the Green Bank Tele¬ 
scope [23]. Before the discovery of the MSP, the system was first identified as a low-mass X-ray 
binary (LMXB) with an accretion disk [24] in 2001. Later, optical and X-ray observations revealed 
the disappearance of the accretion disk, implying a change to the MSP state [25], though the MSP 
was not identified until 2008. However, in 2013 June the radio pulsations disappeared [26], and the 
system was shown to have reformed an accretion disk [27], thus reverting back to the LMXB phase 
and constituting the first time a system has been seen to oscillate between the two states. Further¬ 
more, it has been thought that MSPs are old pulsars that are spun up to rapid rotation frequencies 
through the accretion of material from a companion star [28] in a process sometimes referred to 
as “recycling,” and the observed behavior of the so-called “missing-link” PSR J1023-1-0038 system 
has helped solidify the recycling scenario as the preferred explanantion for the origin rapid rotation 
of MSPs. To our knowledge, this work represents the first time that this intriguing system has 
probed for possible VHE gamma-ray emission. 


2. Observations, Analysis, and Results 


The VERITAS array of four 12 m diameter imaging atmospheric Cherenkov telescopes is lo¬ 
cated at the Ered Eawrence Whipple Observatory (EEWO) in southern Arizona (31° 40'N, 110° 
57' W, 1.3 km a.s.l.) and started full array operations in 2007. The telescope reflectors each consist 
of 345 hexagonal mirror facets, and the cameras comprise 499 photomultiplier tubes giving a total 
field of view of ~3.5°. VERITAS is sensitive to gamma-ray photons in the energy range 0.85 to 
> 30 TeV with a sensitivity to detect a 1% Crab Nebula source in ~25 hr'. It has an energy resolution 
of 15-25%, an angular resolution of 0.1° at 68% containment, and a pointing accuracy error of less 
than 50 arcseconds [29]. 
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2.1 Geminga 



Figure 1: Phase-folded light curves for the Geminga pulsar for both the Fermi-LAT (panels a and b) and 
VERITAS data (panel c). In panel a, the phase-folded light curve of the Ferm/-LAT data with an energy 
threshold of E > 100 MeV is shown, and the region defined as background is indicated by the arrow. Panel b 
shows the data with a higher energy threshold of 10 GeV, and the asymmetric Gaussian fit to P2 is represented 
by the solid black curve. The inset zoom panels for PI and P2 show the respective peaks with energy 
thresholds of both 5 and 10 GeV, with the asymmetric Guassian fit to PI represented by the solid black 
curve. Panel c shows the phase-folded VERITAS data from the location of the Geminga pulsar in units of 
counts per bin, and the average number of counts is indicated by the horizontal dashed line. The vertical 
dashed lines in all three panels represent the phase regions for PI and P2. Eigure from [31]. 


The VERITAS observational campaign on Geminga has resulted in the accumulation of a total 
of 71.6 hr of quality-selected data. Event arrival times are barycentered and phase-folded with 
Tempo2 [30] using an XMM-Newton timing solution (E. Gotthelf 2014, private communication) 
for data obtained before the launch of the Ehnni-EAT. A publically available Fermi-LAI ephemeris 
from a webpage maintained by M. Kerr^ is used for phase-folding all other VERITAS data. The 
phase-folded Geminga pulsar data recorded by the Fermi-LAI shows two emission peaks (PI and 

3www.slac.stanford.edu/~keiTm/fermi_pulsar_timing/ 
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P2), which are used to define phase regions of expected signal for the VHE gamma-ray data taken 
by VERITAS. The peaks PI and P2 are both fitted with asymmetric Gaussian functions above 5 
and lOGeV, respectively, which allows measuring the widths of the pulses at the highest energies 
possible in the Fermi-LPW band. The two expected VHE signal regions are defined a priori as the 
±la regions for PI and P2 ([0.072-0.125] and [0.575-0.617], respectively), and the background 
region is defined as [0.7-1.0] [31]. The phase-folded Femii-LKT and VERITAS light curves are 
shown in Eigure 1 . 

The test for a pulsed signal using the ±la regions defined above reveals no evidence pulsed 
emission from the Geminga pulsar in the VERITAS data. To more generally test for periodicity 
in the phase-folded VERITAS data, the //-Test [32] is employed. The computed H statistic is 
1.8, corresponding to a random chance probability of 0.49, indicating that phase-folded VERITAS 
data is consistent with a random distribution. Upper limits (UEs) on excess counts for the PI 
and P2 phase regions are computed at the 95% confidence level (CE) using the method of [33]. 
These limits are converted to integral flux UEs above 135 GeV of 4.0 x 10^'^cm^^s^^ for PI 
and 1.7 X 10^'^cm^^s^^ for P2. These limits are shown along with a phase-averaged Fermi-LPH 
SED for the Geminga pulsar in Eigure 2. Eor a more complete explanation of this analysis and the 
results, see the recent VERITAS publication on the Geminga pulsar [31]. 


Geminga phase-averaged 
+ SED in 12 Bins [Fermi 5.2 yr] 

- Power-iaw x Exp. Cut-off fit [Likeiihood] 

— — — Power-iaw fit [Likelihood with E > 10 GeV] 
. Power-iaw fit [Extrapoiated above 100 GeV] 

Geminga P1 

□ Upper iimit 95% confidence ievei [VERiTAS] 


Geminga P2 

O Upper iimit 95% confidence ievei [VERiTAS] 

Older Measurments of Geminga 
O Singh et ai, 2009 [PACT] 

▼ Neshpor et ai, 2001 [Crimea] 

V Aharonian et al, 1999 [HEGRA] 

★ Akerlof et al, 1993 [Whipple] 

0 Vishwanath et ai, 1993 [Ootacamund] 
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— ■ — ■ - Crab Nebuia [Albert et ai, 2008] 

— • • • — 1 % Crab Nebuia 

Crab Pulsar [Aliu et al, 2011] 


10 ^ 10 " 10 = 10 = 
Energy [MeV] 



Figure 2: Phase-averaged SED for the Geminga pulsar derived from 5.2 yr of Fermi-LAI data with VHE 
flux UEs (open squares and circles) for the PI and P2 phase regions computed from the VERITAS data. 
The energy thresholds for the two sets of VHE flux UEs are 135 and 550 GeV. Eor a full description of the 
Fermi-LAIT data analysis, see the recent VERITAS collaboration publication on the Geminga pulsar [31]. 
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2.2 PSR J1023+0038 


J1023 State 

H statistic 

2fT pulsed VHE 

3 O' pulsed VHE 



fluxUL(m-2 s-‘) 

fluxUL(m-2 s-‘) 

Radio MSP 

0.28 

6.16 X 10^'^ 

9.73 X lO^** 

Accretion/LMXB 

0.18 

1.12x 10-* 

1.97 X 10^* 


Table 1: H statistics and integral pulsed VHE flux upper limits > 166 GeV for both the radio MSP and 
accrection/LMXB states of PSR 11023+0038. 

The search for pulsed VHE emission 
was split into two parts as VERITAS data 
were taken both before and after the 2013 
June disappearance of the radio pulsar 
(henceforth these two states of the system 
are referred to as the ‘radio MSP’ and ‘ac- 
cretion/EMXB’ states). The VERITAS live 
time accumulated on the target is 18.1 hr for 
the radio MSP state and 8.2 hr for the accre- 
tion/EMXB state. The data span two con¬ 
figurations of the VERITAS array: before 
and after the upgrade of the camera photo- 
mulitplier tubes in the Summer of 2012 [34]. 

Events are phase-folded in Tempo2 using a 
Jodrell Bank radio ephemeris (A. Archibald, 
private communication), though the location 
of phase zero is unknown, precluding the 
possibility of defining signal regions in the 
phaseogram. Therefore only the //-Test is 
used to test for the presence of a periodic 
signal in the data. The //-Test does not re¬ 
veal any evidence for periodicity in the VHE 
phaseograms, which are shown in Eigure 3. 

The H statistics are used to compute 2 and 
3a upper limits via the method of [35] above 
an energy threshold of 166 GeV assuming 
a duty cycle of 10% and Gaussian pulse 
shapes. The computed H statistics and pulsed flux UEs are shown in Table 1. A search for 
steady emission > 300 GeV in the data is also performed, giving significances of —0.8a and 0.2a 
for the radio MSP and accretion/EMXB states, respectively. Upper limits on a steady emission 
component are computed at the 95% CE for both states assuming a photon index of E = 2.5 us¬ 
ing the Rolke method [36]. The computed steady flux 95% CE UEs are 8.1 x lO^^m^^s^^ and 
9.6 X lO^^m^^s^' for the radio MSP and accretion/EMXB states, respectively. These results are 
part of a VERITAS collaboration publication that is currently in preparation. 



Figure 3: Phase-folded light curves of events 
recorded by VERITAS from the direction of PSR 
J1023+0038. One rotation of the pulsar is shown. 
The top panel shows the data collected during the ra¬ 
dio MSP state of the system, while the bottom panel 
shows the data collected during the accretion/EMXB 
state. The solid and dashed lines represent the average 
and error on the average number of counts, respec¬ 
tively. 
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3. Future Plans 

Pulsar science with VERITAS is an ongoing effort and will continue in the coming years. 
There are currently 19 pulsar locations for which VERITAS data exist, including coverage of the 
top ten Fermi-LKV detected pulsars in the Northern Hemisphere when ranked in spin-down power 
divided by distance squared (for details, see Archer, A., these proceedings) [37]. The analyses of 
these archival data sets is currently underway and will represent the first time that any of these 
pulsars have been probed for possible pulsed VHE emission. 

Due to the good coverage of the young pulsars in archival VERITAS data, VERITAS is begin¬ 
ning to target a different promising population of pulsars to probe for VHE emission: millisecond 
pulsars. Because of the rapid rotation of MSPs (spin periods reaching down to ~ 1 ms), these pulsars 
have much smaller light cylinder radii, implying more compact magnetospheres. The compactness 
of their magnetospheres suggests significantly different physical parameters of the emission mech¬ 
anism compared to the young pulsar population. Eurthermore, though less luminous, their HE 
gamma-ray producion efficiencies have been found fo be comparable or greafer fhan fhose of fhe 
young pulsars [2]. Eor fhese reasons, fargefing MSPs in addition fo young pulsars should provide 
complemenfary benefifs for consfraining models affempling fo explain fheir gamma-ray emission. 
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